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SECTION  I 


INTRODUCTION 


During  the  past  three  years,  we  have  been  studying  the  applica¬ 
tions  of  free  electron  lasers  for  Navy  missions  as  well  as  the  physics 
and  technical  issues  of  the  low  voltage  free  electron  lasers.  By  low 
voltage  we  imply  electron  beams  of  kinetic  energy  <  10  MeV.  In  these 
devices  the  conventional  magnetostatic  wiggler  is  replaced  by  an  elec¬ 
tromagnetic  pump  of  wavlength  <  1  mm  to  generate  doppler  shifted  co¬ 
herent  radiation  in  the  visible  portion  of  the  electromagnetic  spec¬ 
trum.  In  our  last  study,  we  found  that  the  extraction  efficiency  and 
available  power  with  an  electromagnetic  pump  depend  critically  on  the 
emittance  of  the  electron  beam.  We  found  further  that  the  losses  in 
the  optical  cavity  of  the  pump  beam  must  be  lO'’1*  per  round  trip  if 
the  overall  system  efficiency  is  to  be  acceptable. 

In  this  study,  we  look  at  several  methods  of  reducing  electromag¬ 
netic  pump  optical  cavity  losses.  The  only  technology  that  seems  to 
show  any  promise  of  being  able  to  achieve  the  required  low  optical 
losses  is  the  use  of  superconducting  cavity.  Even  here  there  are 
physics  and  technology  limitations.  In  general  the  system  works  better 
at  pump  wavelengths  >  1  mm.  In  Section  II,  we  summarize  the  results  of 
the  systems  requirement  analyses  on  cavity  optics  for  the  two  stage 
FEL.  In  Section  III,  we  discuss  in  detail  the  several  methods 
suggested  for  decreasing  cavity  losses.  Section  IV  contains  several 
design  maps  for  choosing  the  optical  parameters  of  the  system.  The 
section  contains  a  summary  of  parameters  for  a  point  design  which  is 
geared  to  the  near  term  experiment  of  UCSB.  Section  V  is  a  summary  of 
the  optical  analysis  reported  with  conclusions  and  recommendations. 
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SECTION  II 


OPTICAL  CAVITY  REQUIREMENTS 

In  this  section,  we  summarize  the  results  obtained  by  us  in  the 
previous  study^.  It  has  been  shown  that  the  constant  parameter  mag¬ 
netostatic  wiggler  can  result  in  an  energy  extraction  of  ~  1/2  N  per 
pass  where  N  is  the  number  of  periods  provided  the  FEL  is  homogeneously 
broadened.  For  the  case  of  inhomogeneous  broadening,  the  extraction 
efficiency  is  even  lower.  To  improve  the  extraction  efficiency  in  a 
single  pass,  it  was  proposed  to  vary  the  parameters  of  the  wiggler  in 
such  a  way  that  the  ponderomot ive  potential  well  in  which  the  electrons 
move  is  decelerated  and,  thereby,  decelerating  the  electrons.  In  a 
magnetostatic  wiggler,  both  the  period,  as  well  as  the  wiggler  magnetic 
vector  potential,  can  be  designed  to  vary  in  a  prescribed  manner. 
However,  with  an  electromagnetic  pump  of  wavelengths  less  than  1  mm, 
first,  the  pump  vector  potential  is  too  small  to  contribute  signifi¬ 
cantly  to  the  resonance  condition;  and  secondly,  it  is  not  possible  to 
change  the  wavelength  of  the  pump.  It  has  been  suggested  that  the 
electron  beam  can  be  made  to  interact  at  varying  interaction  angle  to 
the  pump  radiation  (instead  of  tt  radians)  such  that  the  effective  pump 
wavelength  is  altered.  Such  a  scheme,  though  attractive,  would  severe¬ 
ly  restrict  the  interaction  length.  An  alternate  (and  more  practical) 
scheme  is  to  use  an  axial  accelerating  field  to  keep  the  electrons  at 
resonant  energy  and  to  keep  them  trapped  in  the  ponderomot ive  potential 
well.  For  effective  trapping,  the  electrostatic  coulomb  force  between 
the  electrons  has  to  be  much  smaller  than  the  attractive  ponderomot ive 


potential  force.  In  this  regime  of  operation,  then,  single  particle 
physics  is  applicable. 

Neglecting  the  energy  supplied  by  the  accelerating  field,  the 
extraction  efficiency  is  given  by 


ext 


'ext 


(y  “  1)  mc^  J 


(2.1a) 


eb 


2.68  x  10~42  P3  X  X4  L 

_ UP  ,PP 

(Y  ~  1)  or  C3  J 


(2.1b) 
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where  P  is  the  pump  power  in  watts,  X  and  X  are  the  pump  and  out- 
p  pp  sp  r 

put  wavelengths  in  microns,  L  is  the  interaction  length  in  cm,  J  is 

eb 

the  current  density  in  amp/cm2,  is  the  out-coupling  coefficient  and 

a  is  the  factor  that  determines  beam  quality. 

In  Figure  (1),  the  extraction  efficiency  is  plotted  as  a  function 

of  pump  power  for  various  values  of  the  parameter  a  while  fixing  the 

equivalent  spread  due  to  emittance  at  0.1%.  This  meant  fixing  a2j  v  at 

eb 

~  64  amp/cm2.  The  figure  also  points  out  the  important  gain  that  can 
be  achieved  in  improving  the  beam  emittance.  The  horizontal  hatched 
line  indicates  a  trapped  fraction  of  40%  which  is  about  the  maximum 
that  one  can  achieve  in  a  nonadiabatic  accelerating  field.  This  is 
also  the  cold  beam  limit. 

The  overall  system  efficiency  of  a  low  voltage  free  electron  laser 
that  utilizes  an  electromagnetic  pump  can  be  written  as 


_  _ _  c  acc _ _ _ 

sys  (1-0  I  [  (Y~l)mcz  +  (1-n  )  E  Lj  + 
R  L  c  acc  J 

where  n  is  the  recovery  efficiency  of  the 

the  round  trip  fractional  loss  of  the 


n  e  Li+gp/n 
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electron  beam  energy, 
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Figure  1.  Effect 


efficiency  of  the  pump  laser  system.  In  Figures  2  and  3,  we  have 
plotted  the  system  efficiency  versus  3  for  a  representative  set  or 
parameters  with  as  the  third  variable.  It  is  clear  from  the  figure 

that  for  devices  of  interest  to  Navy  (medium  to  high  powers)  that  the 
round  trip  losses  of  the  pump  beam  must  be  less  than  for  nnv 
reasonable  system  efficiency.  We  are  thus  forced  to  design  reflective 
optics  for  the  pump  cavity  with  >  99.9%  reflectance  per  surface. 

For  output  powers  in  excess  of  tens  of  kilowatts  in  the  visible/- 
near  IR,  the  circulating  pump  power  required  is  in  excess  of  10^ 
watts,  while  the  circulating  power  at  the  short  wavelength  would  be  in 
excess  of  106  watts.  One  is  then  faced  with  the  problem  of  designing 
and  developing  optical  cavities  at  two  disparate  wavelengths  with  high 
to  very  high  reflectivity  at  both  the  wavelengths.  The  optical  design 
could  be  somewhat  simplified  if  means  can  be  found  for  separating  the 
two  wavelengths  so  that  one  can  concentrate  on  the  design  of  high 
reflectivity  optics  at  a  single  wavelength. 


Figure  3. 


System  Efficiency  vs.  Pump  Cavity 

Loss  at  High  Current 


SECTION  III 


LOW  LOSS  CAVITIES 

3 . 1  Room  Temperature  Metal  Reflectors 

The  losses  in  an  electromagnetic  cavity  can  be  classified  as  (i) 
absorption  losses  on  the  cavity  walls,  (ii)  diffraction  losses  in  open 
cavities  and  (iii)  out  coupling  losses.  The  diffraction  losses  can  be 
almost  completely  eliminated  by  going  to  closed  cavities  (save  for 
letting  e-beam  in  and  out).  However,  as  we  shall  see  later,  even  in 
open  resonators,  they  can  be  kept  at  10“4  to  10“5  per  round  trip  and 
thus  would  not  pose  significant  problems.  The  out  coupling  losses  come 
about  because  the  output  radiation  at  the  short  wavelength  has  to  be 
coupled  out  and  the  pump  radiation  will  leakout  at  this  place.  By  the 
proper  choice  of  the  transverse  cavity  mode,  this  can  be  minimized. 
The  most  significant  loss  then  would  be  the  absorption  loss  on  the 
cavity  wall  itself.  In  this  section,  we  shall  assess  the  possible 
methods  of  minimizing  the  absorption  loss  based  on  (a)  material 
property  and  (b)  on  the  geometry  of  the  reflecting  surface. 

With  conventional  mirrors,  one  can  consider  (a)  pure  metallic 
conditions  and  (b)  dielectric  coatings  for  increasing  the  reflectivity. 
For  good  conductors,  the  surface  resistance  Rs  can  be  written  as4 

R  »  l/o5  (3.1) 

s 

where  a  is  the  conductivity  and  6  is  the  skin  depth.  The  skin  depth  d 
is  given  by 

5  =  (2/ojop)1/2  (3.2) 
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where  U)  is  the  frequency  of  the  wave  and  p  is  the  permeability  of  the 
metal.  The  absorption  loss  at  the  surface  of  the  metal  can  then  be 
written  as 

oc  =  R  /Z  (3.3) 

L  so 

where  Z0  is  the  free  space  impedence  (377  ohms).  Table  I  gives  the 
surface  resistance,  skin  depth  and  loss  factor  for  different  metals  at 
500  u  and  1000  p.  The  loss  factor  given  in  the  table  should  be 
multiplied  by  2  to  obtain  6  for  the  cavity  since  two  surfaces  are 
involved.  From  this,  we  find  that  the  best  that  one  can  achieve  in  the 
500  to  1000  p  wavelength  range  is  to  have  a  round  trip  absorption  loss 
of  0.1%.  The  values  calculated  in  Table  I  are  based  on  the  assumption 
of  idealized  surface  conditions.  Surface  contamination  and  imperfec¬ 
tions  on  grain  boundaries  typically  increase  the  absorption  by  a  factor 
of  2. 

The  dielectric  coatings,  while  increasing  reflectivity  signifi¬ 
cantly  in  the  visible,  do  not  offer  the  same  advantage  in  the  far 
infrared.  In  practice  it  is  very  difficult  to  improve  over  the  reflec¬ 
tivity  of  silver^.  The  only  problem  with  silver  seems  to  be  the  oxide 
formation  on  the  surface  (tarnishing).  Workers  in  the  field  prefer 
gold  plated  mirrors  to  silver  for  this  reason.  However,  for  the  appli- 
caton  that  we  have  in  mind,  silver  would  seem  to  be  the  best  choice. 

3 . 2  Reflection  of  an  Electromagnetic  Plane  Wave  by  an  Infinite 
Set  of  Parallel  Wires 

In  this  section  we  discuss  the  possibility  of  decreasing  the 
mirror  losses  by  replacing  it  with  a  set  of  parallel  wires  spaced  less 
than  X/2  apart.  Such  an  array  of  wires  or  tubes  is  frequently  used  in 
the  microwave  region  to  attenuate  waves  passing  through  such  a  struc¬ 
ture.  The  configuration  that  we  have  in  mind  is  shown  in  Figure  4. 
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Mumfordb  has  given  an  eraperical  formula  for  the  ratio  of  incident  to 


transmitted  power  (P^/P^)  for  such  a  geometry. 


where 


(3.4) 


I  B  I 


A  .  .  r  0.83  exp  (2rrr/a)  -j 
a  '  L  exp  (2nr/a)  -1  J  ‘ 


(3.5) 


The  electric  field  vector  is  parallel  to  the  direction  of  the  wires  and 

the  propagation  direction  is  normal  to  the  plane  of  the  wires.  This 

relationship  is  shown  plotted  in  Figure  5  for  several  values  of  the 

ratio  a/2r.  A  word  of  caution  is  in  order  about  equation  (3.5).  It 

does  not  go  to  the  right  limits  when  a*2r .  Thi~  is  to  be  expected 

since  equation  (3.5)  is  only  an  empirical  fit,  but  the  author  claims 

an  accuracy  of  1  db  over  interesting  ranges  of  a/2r  (i.e.  from  3  to 

30).  In  order  to  show  significant  improvements  over  commercially 

available  mirrors  one  would  like  to  have  P.  /Pm  >  10^.  One  would 

inc  Trans  ~ 

also  like  to  have  a/2r  as  large  as  possible  since  the  formula  doesn 1 1 
include  absorption  losses  by  the  wires.  Fabrication  technology  would 
limit  the  smallness  of  the  radius  of  the  wires.  Table  II  gives  the 
wire  diameter  necessary  to  achieve  a  transmission  loss  of  40  db  for 
different  values  of  a/2r.  Whether  this  would  be  possible  remains  to  be 
seen.  In  addition,  the  fiuence  level  that  can  be  withstood  by  the 
wires  would  be  significantly  lower  than  that  for  the  plane  mirror  which 
can  be  cooled  more  easily.  We  are  therefore  forced  to  conclude  that 
this  may  not  be  a  practical  approach  to  making  high  reflectivity 
cavities . 
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TABLE  II 


WIRE  DIAMETERS  CORRESPONDING  TO  PTU<,/Pn.DAkl  =  104 

INC  1  RAN 


Diameter  =  2  r 

A  urn 

a/2r 

Um 

1000 

3 

6.8 

1000 

4 

_  -  | 

3 

1000 

10 

0.44 

3 . 3  Reflection  of  an  Electromagnetic  Plane  Wave  By  an  Infinite 
Set  of  Plates 

It  has  also  been  suggested  that  the  losses  in  a  plane  mirror  can 
be  reduced  by  changing  the  geometry  to  a  set  of  parallel  plates  spaced 
<  A/2  apart  such  that  waves  with  electric  field  vector  polarized 
parallel  to  the  plates  cannot  propagate  inside  the  plates  (evanescent 
wave  only).  This  geometry  is  illustrated  in  Figure  6.  For  the  case  of 
perfect  conductor,  the  set  of  plates  act  as  perfect  plane  mirror.  In 
the  case  of  a  good  conductor,  the  hope  was  that  the  absorption  losses 
would  be  significantly  reduced  from  that  of  the  plane  mirror  because  of 
the  geometrical  reduction  in  the  area  of  the  arrays  normal  to  the 
direction  of  incidence  of  the  electromagnetic  wave.  In  this  section, 
we  shall,  however,  prove  that  this  is  not  the  case.  When  A  <  2a  where 
a  is  the  separation  distance  between  the  plates,  significant  amount  of 
field  is  set-up  between  the  plates  even  though  it  decays  with  distance. 
The  surface  current  set-up  by  this  field  actually  causes  more  power  to 
be  dissipated  than  in  the  case  of  a  plane  mirror. 

The  problem  of  the  field  distribution  when  a  plane  electromagnetic 
wave  is  incident  upon  an  infinite  set  of  staggered,  equally  spaced, 
semi-infinite  metallic  plates  of  zero  thickness  and  perfect  conduc¬ 
tivity  was  solved  in  a  set  of  papers  by  Carlson  and  Heins  in  1947.  The 
boundary  value  problem  for  the  field  was  formulated  mathematically  as 
Weiner-Hopf  integral  equation  and  a  rigorous  solution  to  the  integral 
equation  was  obtained  by  Carlson  and  Heins  by  Fourier  transform 
techniques.  An  expression  for  the  surface  current  density  and  hence 
the  field  was  obtained  in  the  Fourier  transform  domain.  The  actual 


inversion  of  this  transform  to  the  physical  space  was  not  carried  out 


by  the  authors.  Nevertheless,  they  were  able  to  obtain  asymptotic 
expression  for  the  reflected  and  transmitted  fields.  We  are  however 
able  to  use  their  results  and  obtain  closed  form  expression  for  the 
power  dissipated  in  such  an  array  of  plates  when  we  relax  the  infinite 
conductivity  requirement.  The  current  density  J  in  a  good  conductor  is 
confined  to  such  a  small  thickness  just  below  the  surface  of  the  con¬ 
ductor  that  it  is  equivalent  to  an  effective  surface  current  K 

ef  f 

Keff  =  r  j  dC  =  *  X  *11  <3-6) 

where  £  is  the  distance  into  the  conductor  and  is  the  parallel  com¬ 
ponent  of  the  magnetic  field1*.  Comparison  of  equation  (3.6)  with  a 
similar  one  for  the  perfect  conductor  shows  that  a  good  conductor 
behaves  effectively  like  a  perfect  conductor,  with  the  idealized 
surface  current  replaced  by  an  equivalent  surface  current  which  is 
actually  distributed  throughout  a  very  small,  but  finite,  thickness  at 
the  surface.  The  power  loss  per  unit  area  can  be  written  in  terms  of 
the  effective  surface  current: 


dP 


loss 


dA 


2 06  1  Ke  f f 1 


(3.7) 


where  a  is  the  conductivity  of  the  metal  and  6  is  the  skin  depth. 
Thus,  if  we  have  solved  for  the  field  and  currents  for  the  idealized 
problem  of  infinite  conductivity,  we  can  calculate  approximately  the 
absorption  losses  for  the  non-ideal  conductor. 

In  what  follows,  we  shall  be  concerned  only  with  the  case  of 
polarization  such  that  the  electric  field  vector  is  parallel  to  the 
plates.  In  the  case  when  the  H-field  is  parallel  to  the  plates,  the 
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fundamental  mode  can  be  sustained  for  propagation  between  the  plates 
and  hence  is  not  of  interest  for  the  application  that  we  are  discussing 
here.  When  the  E-field  is  parallel  to  the  plates,  all  the  field  compo¬ 
nents  can  be  derived  from  <{>  (x,z)  =  E  (x,z). 

y 

Carlson  and  Heins7  define  a  quantity  I  (z)  which  is  proportional 

YTl 

to  the  surface  current  density  K(z)  on  the  mth  plate.  We  are 
interested  in  the  case  when  the  angle  ct  *  tt/2  and  the  angle  of 
incidence  0=0.  From  Reference  7,  we  have 


J  (w)  = 


2  i  K,  (k  cos  6) 
K_  (w)  (w-k  cos  6) 


(3.8) 


where  J(w)  is  the  Fourier  transform  of  I  (z)  and  w  is  the  transform 

o 

variable.  The  functions  K+  (w)  and  K_  (w)  are  given  bv 


K+(  w) 


(a2  +  b2)(w-q9)  exp  1  x  (W)I  Gi 
2a  G3  (w  +  <)  exp  (iaw/ir) 


(3.9) 


„  ,  N  _  Gtl  (w-<)  exp  [  -i aw/ tt  +  x  (w)] 
K_(w)  "  - (w-op  c5  G6 


(3.10) 


where 


II,  [a  _  i\|>  ]exp  { [kp_wb+wai ) /2irn] +i ( TT/2-a) 
n=l  n  n J 


5*  [A^  +  i^]  exp  {  [  (kp-wb+wai)/2TTn]-i(Tr/2-a)} 


’3  =  J2  l  7  1  -  (  h  )2  ‘  ]  (  7  }  exp 


Gu  =  n  [  V 1  -  (— ^  ]  -  exp  (-iaw/nn) 
4  n=2  L  tt n  TTn  J  TT  r 


Gr  =  H*  [a  -  ity  ]  exp  { 1  (kp-wb+wai)/2im]  +  i(ir/2-a)  } 
5  n  =  -°°  L  n  nJ  r  1 


g6  ~  [^n  +  ]  exp  { [  (kP“wb*-wai)./27Tn]  — i  ( tt / 2— ex ) 

Oj  =  k  cos  6,  Go  =  k  cos  (2a-0), 


/.  9  , ,  kp  ,ak  x2  , ,  kp  x  wa 

sin2  a  (1  -  tt—)  -  (tt—)  ,  ~  cos  a  (1  -  7 5—)  +  - 

2  Tin  9irrr*  Tr*  7lTr> 


'2irn  5 


27Tn 


b  =  a  cot;  a,  p  =  b  cos  6  +  a  sin  0 


X(w)  = 


-law 


[  («  ”  T  )  cot  a  “  In  2  sin  a] 


Since  the  surface  current  density  on  the  plates  is  a  function  of 
distance  into  the  plates,  we  are  really  interested  in 


/"  I K( z ) I  2  dz  =  f  2  f°  I1  (z)l2  dz. 
o  o 

The  quantity  I^(z)  is  defined  to  be  identically  zero  for  z  <  0. 
Therefore,  we  can  use  Parseval’s  formula  and  write 


/  |l  (z)|2  dz  =  /  |l  ( z ) I ^  dz  =  /  )j(w)j2  dw 

J  o  J  o  J 


For  the  special  case  of  a  =  tt/2  and  0=0,  we  find  that 
|  K+  (k  )  |  ^  =  k2a2  and 


[!_(  ii  )2  .  (  )‘] 

|K.(»)(w-k)|2  -  -  9,  *”  "" 

n=l 


wa  v2j 
wa  ^  2 1  2 


f  l-(  *1-  )2  +  (  )0‘ 

1  V  2nTt  V  2mT  '  J 


We  therefore  have, 


/  |j(w)|2  dw  =  -  4  k2  a 


r°°  Jr  n  [1  “  ( p/n) 2  +2(C/n)2] 

i  ^  n-1  [1  -  (2"p/n)z  +  (2C7^FT 


where 


^  =  wa/2ir  and  p  -  ka/2ir  »  a/X  . 


} 


esc  a 
2itn 


.11) 


.12) 


.13) 


-19- 


The  integral  can  be  evaluated  by  the  method  of  residues.  The  poles  of 
the  integrand  are  located  at 


=  ±  i  £  =  -  i  Vn 2/4  - 


(p  <  1/2).  (3.34) 


Evaluating  the  residues,  we  have, 


/  | J(w) I  dw  =  4k2  a  £ 


JI  (l-m2/4n2) 
n=l 


m-1  JL  (l-m2/n^)^  /m2/4-p2 


(3.13) 


The  infinite  products  can  be  written  in  terms  of  simple  trigonometric 
functions  and  the  final  result  is, 


/  |J(w)|2  dw  =  k2  a  l  (2p+l)  2  [ ( 2p+l ) 2  -  4p2]~1/2(3.16) 

00  p  =  0 


The  power  loss  per  sheet  per  unit  distance  in  the  y-direction  is  there¬ 
fore  given  by 

dP  00 

— uSheer  *  •  —  I H  I  2  l  ( 2p+l )  ~2  [  (  2p+l ) 2  -  4p2]~1/2(3.17) 

dy  2o6  tt  r> 

p=0 

In  the  case  of  a  plane  mirror,  the  effective  surface  current  is 


|K  -J  =  In  x  Hi .  (3.18) 

ef  f 

The  power  loss  in  plane  mirror  per  unit  distance  in  the  y-direction  and 
distance  a  in  the  n  direction  is 


d? 
dy  i 

plane  mirror 


(3.19) 
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Comparing  Equation  (3.17)  and  (3.19),  we  can  find  the  ratio  of  power 
loss  with  a  stack  of  sheet  to  that  of  the  plane  mirror.  If  this  ratio 
is  described  by  R,  we  obtain, 


R  =  y  ( 2  p  + 1 )  2  [  (  2p+ 1 )  2  -  4u2]~1/2  (3.20) 

11  P=0 

In  Figure  7  we  have  plotted  R  vs  \i  (*  a/ A).  When  p+0,  we  find  that  R 
tends  to  ~  21.43  while  R*00  as  p-K).5.  The  fact  that  R  is  not  equal  to  1 
at  p=0  is  not  surprising  since  the  theory  was  developed  for  sheets  of 
zero  thickness  and  we  cannot,  therefore,  expect  to  recover  the  limit  of 
plane  mirror  when  the  distance  between  the  sheets  approaches  zero.  In 
fact,  we  do  not  expect  the  results  to  be  strictly  valid  when  the 
distance  between  the  plates  is  less  than  or  comparable  to  the  skin 
depth.  The  conclusion  that  we  can  obtain  from  Figure  7  is  that  the 
idea  of  using  a  stacked  set  of  plates  instead  of  a  plane  mirror  will 
not  result  in  reducing  absorption  losses;  on  the  contrary,  they  will 
increase  the  total  resistance  losses. 

3.4  C jperconduc ting  Cavities 

3.4.1  Introduc  t ion 

From  the  two  preceding  subsections,  it  is  clear  that  the  losses  on 
the  cavity  walls  cannot  be  decreased  significantly  (if  at  all)  by 
purely  geometrical  means.  The  only  other  alternative  to  consider  is 
the  use  of  superconducting  optics  to  reduce  the  losses  significantly. 
There  are  two  fundamental  limits  on  the  operation  of  the  superconduct¬ 
ing  cavities.  These  are  (i)  the  wavelength  limit  and  (ii)  the  surface 
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magnetic  flux  density  limit.  The  experimental  data  on  the  performance 
of  the  superconducting  cavities  lies  mostly  in  the  wavelength  of  a  few 
centimeters  while  we  are  interested  in  millimeter  and  submillimeter 
wavelength  region.  Thus,  caution  must  be  exercised  in  using  the 
extrapolated  values  of  the  surface  resistance  at  these  wavelengths. 

In  a  superconductor  the  lowest  excited  state  for  an  electron  is 
separated  by  an  energy  gap  A  from  the  ground  state.  The  existence  of 
this  gap  in  the  excitation  spectrum  has  been  confirmed  by  a  variety  of 
experiments:  electronic  heat  capacity,  thermal  conduc t ivity ,  ultra¬ 
sonic  attenuation,  far  infrared  and  microwave  absorption,  and 
tunneling.  It  is  to  be  noted  that  the  excitation  of  electrons  across 
the  gap  by  photons  requires  a  minimum  energy  of  2A,  which  ic  consistent 
with  the  description  of  the  superconducting  state  in  terms  of  Cooper 
pairs  as  described  in  Reference  9.  The  energy  gap  decreases  monotoni- 
calLy  with  increasing  temperature  having  a  value  Aq  »  1.75  k  Tc  at 
0°K  and  vanishing  at  the  critical  temperature  Tc  (where  k  is  the 
Boltzmann  constant). 

Thus  for  frequencies  at  which  the  energy  of  the  electrons  hw  is 
smaller  than  2A ,  superconductivity  can  be  achieved.  This  accounts  for 
the  fact  that  at  D.C.  and  low  frequencies  superconductivity  can  be 
achieved  at  temperatures  very  little  below  Tc .  To  operate  at  higher 
frequencies  the  gap  A  must  be  increased. 

The  limiting  frequency  f  at  0°K  for  Niobium  is  717  GHz  (Tc  * 
9.25°K)  while  for  Niobium  Tin  it  is  1395  GHz  (Tc  *  18°K). 

At  temperatures  above  0°K  the  limiting  frequency  decreases  to  zero 
monotonical ly .  No  simple  relationship  for  the  decrease  seems  to  exist. 
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It  is  probably  reasonable  to  assume  as  a  first  estimate  that  one  should 

divide  the  values  of  (f)  by  a  factor  of  2  for  operating  at  a  temper- 

U  K 


ature  above  0°K  and  another  factor  of  2  to  provide  a  margin.  This  then 
would  give  the  limits  as 

f  53  180  GHz  or  X  =  1667  micrometer  for  Niobium 
f  *  349  GHz  or  X  =  860  micrometer  for  Niobium  Tin 
If  the  total  dividing  factor  of  4  is  too  large  it  may  be  possible  to 
reach  somewhat  shorter  wavelength. 

The  surface  resistance  of  a  superconductor  has  been  theoretically 
derived  to  be  equal  to1^ 


rbcs 


a 


1.7 


T 


-A(T)/kT 


where  the  subscript  BCS  refers  to  the  Bardeen,  Cooper  and  Schriefer 
theory.  By  fitting  the  experimental  data  to  the  theory,  one  finds  that 
the  actual  relationship  is 

T 

R  =  2.4  x  10-25  x  — e_1*86  T~  OHMS  (3.21) 
nub  I 

A  complication  arises  due  to  the  fact  that  there  is  a  small  residual 

surface  resistance  R  which  must  be  added  to  R  and  presents  a 

RE  b  BL  b 

lower  limit  which  is  approached  asymptotically  as  T  approaches  0°K. 
It  appears  that  this  residual  resistance  is  very  much  a  function  of 
surface  preparation  and  condition.  The  residual  resistance  may  also  be 
a  function  of  frequency  especially  when  operating  close  to  the  wave¬ 
length  limit  discussed  above. 
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In  order  to  obtain  a  reasonable  extrapolation  to  the  wavelengths 


of  interest  based  on  the  experimental  data  at  10  GHz  “  30  mm  in 

References  10,  11  and  12  we  have  used  the  w**7  scaling  and  these  are 

shown  in  Table  III. 

It  is  seen  in  Table  III  that  Niobium  Tin  at  4.2°K  can  reduce  the 
loss  factor  by  2  to  3  orders  of  magnitude  and  Niobium  at  1.5°K  by  about 
another  order  of  magnitude  as  compared  with  copper  at  room  temperature. 
The  losses  will  in  fact  become  so  low  that  they  will  become  negligible 
compared  with  the  coupling,  stray  and  diffraction  losses.  Note  also 
that  low  losses  make  lower  power  demands  on  the  cryogenic  equipment  and 
ease  thermal  conductivity  oriented  features.  It  is  important  to  recog¬ 
nize  that  only  if  these  other  losses  can  also  be  kept  low  does  it  make 
sense  to  go  to  superconductors. 

In  addition  to  the  wavelength  limit,  a  limitation  on  the  rf  field 
strength  of  superconducting  resonators  exists  and  is  given  by  the 
thermodynamically  critical  flux  density  Bc .  For  type  I  superconduc¬ 
tors  such  as  In.Sn  and  Pb  with  a  low  Bc ,  this  limit  is  really 
observed.  In  type  II  superconductors  external  magnetic  fields  pene¬ 
trate  into  the  superconductor  at  the  lower  critical  flux  density 
Bci»  leading  to  hysteresis  losses.  According  to  more  recent  inves¬ 
tigations,  however,  the  value  Bc^  is  not  important  at  frequencies 
in  the  GHz  range,  so  even  with  type  II  superconductors  the  rf  field 
strength  should  be  limited  by  Bc .  In  general,  however,  the  measured 
rf  field  strengths  are  still  below  this  limit. 

Table  3  of  Reference  10  shows  a  decrease  in  critical  RF  flux 
densities  for  Niobium  cavities  from  1000-1500  Gauss  at  8.6  -  10.5  GHz 
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to  about  800  Gauss  at  3.7  GHz  650-850  Gauss  at  2.8  GHz  and  350  Gauss  at 
1.3  GHz.  For  one  of  the  data  points  the  flux  density  observed  exceeded 
the  lower  critical  flux  density  Bc  ^ .  This  confirms  the  above  claim 
that  the  value  %ci  is  not  necessarily  limiting  at  frequencies  in 
the  GHz  range.  It  remains  to  be  seen  if  this  trend  continues  to 
increase  the  critical  flux  for  the  higher  frequencies  of  interst  which 
would  be  useful  in  our  application.  We  may  take  the  maximum  RF  field 
on  the  surface  of  the  superconductor  to  be  ~  1000  gauss  which  corres¬ 
pond  to  an  intensity  of  ~  10®  w/cm2  on  the  cavity  walls.  If  the  free 
electron  laser  requires  an  intensity  of  ~1010  w/cm2  in  the  interaction 
region,  the  beam  has  to  be  expanded  by  a  factor  of  100  or  greater  in 
area  before  it  hits  the  walls  of  the  superconducting  cavity. 

The  major  problem  with  the  operation  of  the  superconducting  cavity 
is  the  cryogenic  power  requirement.  Typically,  it  takes  1  to  2  kilo¬ 
watts  of  power  to  remove  1  watt  of  heat  at  4°K.  From  the  point  of  view 
of  the  overall  system  efficiency  it  then  becomes  clear  that  nothing 
would  be  gained  by  resorting  to  superconducting  cavity  if  the  losses 
are  not  decreased  by  more  than  three  orders  of  magnitude  compared  with 
room  temperature  cavity.  In  addition  to  heat  addition  to  the  walls  of 
the  superconductor  by  the  pump  radiation,  one  also  has  to  minimize 
radiational  heating  of  the  mirror  by  the  ambient  around  the  mirror. 
One  would  have  to  cool  the  entire  FEL  system  by  successive  stages  of 
cryogenic  cooling.  Our  conclusion  is  thus  that  the  superconducting 
cavity  might  work  in  a  practical  system  if  (i)  the  pump  wavelength  is 
greater  than  a  few  millimeters  and  (ii)  the  operating  temperature  is 
below  ~  2°K.  Clearly,  for  the  proposed  UCSB  experiment  one  would  not 
want  to  complicate  matters  by  using  superconducting  cavity. 
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SECTION  IV 


CAVITY  DESIGN  CONSIDERATIONS 

From  the  preceding  sections,  we  can  conclude  that  for  near  term 
experiments  the  optical  cavity  in  the  far  infrared  is  best  fabricated 
using  silver  mirrors.  For  the  demonstration  of  the  concept  of  two- 
stage  FEL,  it  is  not  necessary  to  have  very  high  reflectivity.  In  this 
section,  we  shall  map  out  the  optical  cavity  parameters  on  a  consistent 
set  of  axes.  A  number  of  assumptions  will  be  made  in  deriving  various 
formulae  and  these  will  be  stated  without  always  a  complete 
justification. 

We  shall  first  assume  that  the  first  stage  wiggler  is  made  out  of 
SmC05  permanent  magnet  as  shown  in  Figure  8.  For  a  magnet  bar  height 
g  =  3  Aw/8  and  a  fill  factor  of  unity,  i.e.  no  gaps  between  neighbor¬ 
ing  magnets,  the  normalized  wiggler  vector  potential  on  axis  can  be 
written13 

a  -  1.07  x  10“4  B  X  exp  [  -  2tt  h/X  ]  (4.1) 

w  r  w  w 

where  Br  is  the  remnant  field  in  gauss  and  other  quantities  are 
defined  in  Figure  8. 

The  photon  beam  width  at  the  entrance  to  the  wiggler  is  denoted  by 
2w  and  for  reasonable  clearance  the  magnet  separation  2h  can  be  taken 
to  be  at  least  a  factor  of  a  larger  than  2w,  The  Rayleigh  range  ZR 
is  related  to  the  photon  beam  waist  radius  wQ  by 

Z  ■  IT  w2  /\  (4-2) 

K  O  , 


JH4 
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The  length  of  the  wiggler  is  taken  to  be  2  q  times  the  Rayleigh 

range.  Thus  the  beam  radius  at  the  entrance  to  the  wiggler  can  he 
written  as 

w  =  w  [  1  +  (  L  / 2Z_) 2 ]  1/2 
o  w  R 

=  [  XZD  (  1  +  q2)/it]  1/2  (4.3) 

Replacing  h  in  equation  (4.1)  by  aw  and  substituting  (4.3)  we  get, 

a  =  1.07  x  10~4  B  X  exp  [  -2TTa  f  A  Z  (1  +  q 2 ) /tt  ]  1 /2/X  ]  (4.4) 

w  r  w  R  wJ 

The  output  wavelength  A  is  related  to  Aw  and  aw  through 

X  =  X  (1  +  a2)/2y2  (4.5) 

w  w 

The  two  convenient  set  of  independent  parameters  are  Aw  and  Zp  and 
we  shall  plot  the  different  physical  quantities  in  maps  with  \w  and 
Zr  as  axes.  Figures  (9)  and  (10)  show  a  plot  of  constant  a^^  lines. 
In  plotting  the  figure,  we  had  to  assume  a  value  for  Y  the  electron 
beam  energy.  The  maximum  value  of  y  that  can  be  obtained  for  the  UCSB 
experiment  is  12  and  we  have  taken  that  to  be  a  fixed  value.  Changing 
Y  would  result  in  similar  plots. 


The  small  signal  gain  per  pass  is  given  by  14 


80LW 

c  2 

r2  n 

1  o  e 

B2  X2 
w 

X3/2  Lw 

(4.6) 

7T 

/  en  mc^ 

(1  +  a*)l/z  Xi/Z  (Au/u)^ 

w  w 

where  r0 

is 

the  classical 

electron 

radius,  ne  is  the 

electron 

density  in  the 

beam  and  (Au/u)  is  the 

fractional  line  width 

l  and  is 

taken  to 

be  ~ 

1/2N  (N  being 

the  number 

of  wiggler  periods). 

If  we 

assume  that  the  photon  beam  waist  radius  and  electron  beam  radius  are 
matched,  equation  (4.6)  can  be  rewritten  as 
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8.70  x  10~3  ( 


(4.7) 


S0Lw 


3w  ZK  q 


(  )/  amp 


where  I  is  the  total  electron  beam  current.  The  saturation  flux  for  a 
constant  parameter  wiggler  is  given  by 


4> 


sat 


2.67  x  107  (  -  )4 

q 


A2 

w 


a* 

w  R 


(4.8) 


For  a  given  set  of  Aw  and  ZR,  a^,  can  be  determined  from  eq. 

(4.4),  Using  eqs.  (4.7)  and  (4.8),  we  can  then  find  the  small  signal 

gain  and  the  saturation  flux.  Figures  (11)  to  16)  show  contours  of 

constant  g  L__,  <{>  and  output  wavelength  plots  for  q  =  1  and  q  =  2. 
o  W  s  a  t 

In  all  the  cases  the  magnet  gap  to  optical  beam  diameter  at  entrance 

and  exit  to  the  wiggler  was  taken  to  be  equal  to  2.  For  the  UCSB 
experiment,  the  same  electron  beam  is  planned  to  be  used  for  the  second 
stage  also.  Thus,  if  we  fix  Y  at  12  and  pick  the  first  stage 
parameters,  the  second  stage  output  wavelength  is  predetermined.  These 

are  also  indicated  in  Figure  (16).  Now  we  are  in  a  position  to  specify 

the  parameters  for  the  proposed  UCSB  two  stage  experiment.  The  Van  de 
Graaff  generator  is  expected  to  provide  10  amperes  current.  The 
bandwidth  requirement  on  the  pump  wave  to  be  a  coherent  wave  for  the 
second  stage  is  that  AX/ A  ~10~4.  The  spontaneous  emission  from  the 
magnetic  wiggler  is  given  by15 


=  2.94  x  10  ®  .  hv  .  i  .  N  watts/ev  of  bandwidth 
dQ  ev  ma 


2 

/ma/mi 1 1 irad ian  /cm 


(4.9) 
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Figure  11.  Map  of  Constant  Small  Signal  Gain 
Contours  for  Higgler  Length  Equal 
to  Twice  the  Rayleigh  Range 


Figure  12.  Map  of  Constant  Small  Signal  Gain 
Contours  for  Wiggler  Length  Equal 
to  Four  Times  the  Rayleigh  Range 


For  a  pump  wavelength  of  230v  and  N  =  5,  the  spontaneous  intensity 


within  the  bandwidth  of  interest  is  found  to  be  ~  6.7  x  i0~® 
watts/cm^.  One  would  need  ~  35  exponentiations  to  get  to  saturation 
levels  of  10®  to  10^  w/cm^. 

In  Figure  (17)  we  have  plotted,  once  again,  lines  of  constant 
gain,  cavity  flux,  output  purnp  wavelength  and  in  plane.  For 

y  =  12,  the  pump  wavelength  is  limited  to  230  to  460p  if  we  stipulate 
that  the  output  wavelength  is  in  the  visible  (4000  to  8000  A° ) .  One 
would  also  need  a  small  signal  gain  of  at  least  5%  per  pass  (requiring 
~  700  passes  for  saturation).  This  is  important  as  it  directly  impacts 
on  the  pulse  length  necessary  for  a  successful  two  stage  operation. 
From  the  point  of  view  of  wiggler  fabrication,  one  would  want  to  keep 
aw  1.  To  achieve  reasonable  second  stage  gain,  we  should  have  the 
cavity  flux  greater  than  ~  10®  w/cm^.  The  shaded  area  in  Figure 
(17)  indicates  the  available  parameter  space  for  operation.  The  actual 
point  that  one  would  pick  in  this  region  depends  on  the  pulse  length  of 
the  electron  beam  available.  If  we  take  the  cavity  length  to  be  ~  3 
meters,  one  round  trip  time  is  20  nsec.  Thus  for  700  passes,  the  pulse 
length  necessary  is  14  Msec.  To  reduce  the  pulse  length,  we  have  to 
pick  an  operating  point  with  higher  small  signal  gain  at  the  expense  of 
decreasing  saturation  flux.  This  would  reduce  the  second  stage  gain 
and  so  we  cannot  arbitrarily  increase  the  small  signal  gain.  We  can, 
however,  locate  the  first  stage  interaction  region  away  from  the  beam 
waist  and  hope  to  operate  the  second  stage  interaction  region  at  the 
beam  waist.  Usually,  the  advantage  gained  by  such  tricks  is  not  over¬ 
whelming.  If  we  choose  g0L^  to  be  10%  per  pass  (for  10  ampere 
current),  then  for  a  3  meter  cavity,  we  would  need  ~  7  psec  long  pulse 
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(neglecting  optical  losses).  A  2  meter  cavity  would  reduce  this  to  ~  5 


Usee.  Figures  (18)  to  (20)  show  similar  plots  as  Figure  (17)  but  for 
different  values  of  q  and  y.  Using  these  figures,  we  can  pick  a  con¬ 
ceptual  design  point  for  the  operation  of  the  two  stage  FEL.  It  must 
he  borne  in  mind  that  the  boxed  region's  boundary  is  somewhat  fuzzy  and 
arbitrary.  The  requirement  on  the  output  wavelength  can  be  consider¬ 
ably  relaxed  if  one  uses  two  separate  electron  beam  sources  for  the 

first  and  second  stage.  The  limit  on  the  minimum  gain/amp  can  be 
relaxed  if  the  current  can  be  increased  to  say  100  mperes.  Comparing 
Figures  (17)  and  (18),  we  see  that  q  =  2  gives  a  bigger  operating 

window  than  q  =  1 .  One  does  not  want  to  arbitrarily  increase  q  and 

operate  with  a  smaller  Rayleigh  range  because  (a)  the  photon  beam  size 
at  the  waist  would  approach  the  pump  wavelength  and  (b)  the  formulas  we 
have  used  here  assumed  sensibily  constant  flux  over  the  length  of  the 
wiggler  and  therefore  would  not  be  applicable.  For  the  conceptual 

design  of  the  two  stage  FEL,  we  would  therefore  choose  the  nominal 
parameter  values  given  in  Table  IV  for  the  UCSB  experiment. 

For  a  10  Msec  long  pulse,  the  fluence  at  the  beam  waist  is 

~  6  KJ/cm2.  An  area  expansion  of  ~  100  would  bring  the  mirror  fluence 

to  60  J/cm2  and  energy  absorbed  with  0.2%  absorption  would  be 

~  120  mj/  cm2,  which  should  be  quite  safe  interns  of  mirror  damage.  If 
Lc  is  the  length  of  the  optical  cavity  and  Rj  and  R2  are  the 


mirror  radii,  we  can  define  quantities  gj  and  g2  by 

gl>2  =  1  (^.10) 

For  stability16,  we  require  that  0<  gjg2  <  1.  We  shall  consider  the 
symmetric  case  with  gj  *  g2  *  g*  Since  we  want  the  beam  to  expand 

rapidly,  we  shall  choose  a  design  where  the  cavity  is  nearly  concentric 
(L  *  2R) . 


TABLE  IV 


WIGGLER  PARAMETERS  FOR  THE  TWO-STAGE  FEE  EXPERIMENT 


Y 

= 

12 

I 

= 

1 0  amp 

X 

— 

4  cm 

w 

ZR 

* 

6  cm 

Lw 

> 

24  cm 

aw 

= 

0.855 

8oLW 

= 

9.54%  per 

pass 

*s 

= 

5.8  x  LO8 

w/cm^ 

w 

= 

2 . 14  mm 

0 

X 

= 

240  u 

P 

X 

5= 

419  nm 

s 

Number 

of 

periods  - 

6 
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Let  L/R  -2-e.  Then 

g  “  e  -  1  . 

The  radius  of  the  beam  at  the  waist  is  given  by  16 


w 

o 


1/2 


XL 

=  (  — -  ) 
V  2tt 


1/2 


Lju. 

[  8! 


1/4 

IZ-LL  -..&ug.2lrrz 

+  §2  -  2^1  82 5 

1/4 

) 


Equation  (4.11)  can  be  rewritten  as 


71  W 


XZD  L 
R  c 


(  ) 

ZR  2-c 


1/2 


Since  Z„  is  defined  to  be  tt  w2/X,  we  have, 
R  o 


q  [  e/ (2-e) ] 
c 


1/2 


=  2 


(4.11) 


(4.12) 


(4.13) 


where  q  »  L  /Z_. 
c  c  R 

The  spot  size  on  the  mirrors  is  given  by 


XL 


1/2 


”1  -  <  ~  > 


1/4 


81  (1  _  81  g2> 


This  can  be  cast  in  the  form, 

q  [  - ~ - —  1  =  w2/w2 

\  1  e  (2-e)  J  1  o 


From  equation  (4.13)  and  (4.15),  we  obtain, 


e  =  2  w2/w2 
o  1 


O  ,  2-e  1/2 

<L  "  2  (  T-  > 


(4.14) 


(4.15) 


(4.16) 

(4.17) 


Knowing  qc  and  e,  we  can  obtain  the  cavity  length  and  mirror  radius 
of  curvature.  The  diameter  of  the  mirror  can  be  calculated  by  fixing 
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the  diffraction  toss. 


Lii7  has  numerically  calculated  the  diffraction 
losses  for  different  |g|  values,  and  Fresnel  numbers.  In  general,  if 
Igl  <  0.99,  one  would  need  Fresnel  number  of  15  or  larger  to  keep  the 
diffraction  losses  per  pass  to  <  10~4  .  In  Table  V  two  different  cavitv 
designs  are  given  for  an  area  expansion  factor  of  100  and  200  for  a 
Fresnel  number  of  15  (=  a^/ALc,  where  a  is  the  radius  of  the 
mirror) . 

The  choice  between  the  first  and  second  design  would  depend  on 
other  factors  such  as  space  to  bend  the  electron  beam  into  and  out  of 
the  wiggler,  electron  beam  pulse  length,  etc. 
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SECTION  V 


SUMMARY 

We  have'  invest igated  several  methods  of  reducing  the  pump  cavity 
losses  in  a  two  stage  free  electron  laser.  We  have  found  that  with 

metal  optics,  we  can  reduce  the  round  trip  cavity  absorption  losses  to 
0.2%  with  some  difficulty.  The  cavity  wall  losses  can  be  signifi¬ 
cantly  decreased  by  resorting  to  superconducting  optics.  The  limita¬ 
tions  of  this  system  is  that  the  operating  wavelength  has  to  be  a  few 
millimeters  and  the  operating  temperature  has  to  be  <  2°K  to  keep  the 
refrigeration  power  at  a  reasonable  level.  One  has  also  to  pay  the 

price  of  the  complexity  of  the  system.  For  near  term  UCSB  experiment, 
it  was  felt  that  room  temperature  silver  coated  optics  would  be  the 
best  choice.  To  expand  the  optical  beam  to  low  flux  levels  that  the 
mirrors  can  handle,  it  was  necessary  to  use  near  concentric  resonator 
geometry.  For  such  a  configuration,  the  diffraction  losses  could  be 
maintained  <10~4  by  choosing  Fresnel  number  greater  than  or  equal  to 
15.  <  ‘Oool 

A  permanent  magnet  Samarium  Cobalt  wiggler  configuration  was 
chosen  for  the  opt imizat Lon  of  the  wiggler  design  procedure.  We  found 
that  all  the  key  parameters  for  a  successful  design  of  the  two  stage 
FEL  could  be  mapped  in  the  convenient  parameter  space  of  wiggler  period 
and  Rayleigh  range  of  the  pump  photon  beam.  By  overlaying  the  several 
parameters  such  as  small  signal  gain  per  pass,  saturation  flux,  output 
wavelength  etc.  in  the  same  parameter  space,  we  were  able  to  identify  a 
region  for  the  operation  of  the  first  stage  of  the  two  stage  FEL.  This 
enabled  us  to  pick  a  reasonable  set  of  FEL  parameters  for  the  first 
stage.  Knowing  the  bearr  waist  radius  and  saturation  flux,  we  could 
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calculate  the  spot  size  on  the  mirror  and  expansion  needed  to  handle 
the  optical  flux.  This  in  turn  enabled  us  to  obtain  the  cavity  length 
and  mirror  radii  of  curvature.  Using  a  Fresnel  number  of  15,  we 
finally  obtained  the  mirror  size. 

Detailed  analytical  calculations  have  to  be  carried  out  to  calcu¬ 
late  saturation  flux  and  other  parameters  correctly.  This  in  turn 
would  impact  on  the  optical  cavity  design  slightly.  The  optimization 
procedure  discussed  gives  a  range  of  design  parameters  (albeit  narrow) 
for  a  given  output  wavelength.  Further  optimization  may  be  carried  out 
depending  on  the  applications  requirements,  etc. 
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